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ABSTRACT 



Models that reproduce the observed high-velocity clouds (HVCs) also predict 
clouds at lower radial velocities that may easily be confused with Galactic disk 
{\z\ < 1 kpc) gas. We describe the first search for these low- velocity halo clouds 
(LVHCs) using IRAS data and the initial data from the Galactic Arecibo L-band 
Feed Array survey in Hi (GALFA-Hl). The technique is based upon the expecta- 
tion that such clouds should, like HVCs, have very limited infrared thermal dust 
emission as compared to their Hi column density. We describe our 'displacement- 
map' technique for robustly determining the dust-to-gas ratio of clouds and the 
associated errors that takes into account the significant scatter in the infrared 
flux from the Galactic disk gas. We find that there exist lower-velocity clouds 
that have extremely low dust-to-gas ratios, consistent with being in the Galactic 
halo — candidate LVHCs. We also confirm the lack of dust in many HVCs with 
the notable exception of complex M, which we consider to be the first detection 
of dust in HVCs. We do not confirm the previously reported detection of dust in 
complex C. In addition, we find that most Intermediate- and Low- Velocity clouds 
that are part of the Galactic disk have a higher 60/im/lOO /im flux ratio than 
is typically seen in Galactic Hi, which is consistent with a previously proposed 
picture in which fast-moving Galactic clouds have smaller, hotter dust grains. 

Subject headings: ISM: kinematics and dynamics, ISM: clouds. Galaxy: halo, 
galaxies: formation 



1. Introduction 



High- Velocity Clouds (HVCs) are clouds of neutral Hydrogen (Hi) in or around the 



Galaxy that are inconsistent with Galactic rotation. Since the discovery of HVCs by MuUer 
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et al. (1963) many theories have been proposed for their origin (e.g. Oort 1966), including 



Galactic fountain models (e.g. Bregman 1980), material stripped from accreting satellite 



galaxies (e.g. Mayer et al. [2006 ) and formation in the cooling Galactic baryonic halo 



e.g. Mailer & Bullock 2004). These origins can be constrained by the observed distribution 



of HVC metallicities, distances, morphologies, fluxes, and velocities. Indeed, many of the 
originally proposed origins have been disregarded as they have failed to meet some of these 
observational criteria. While many of these models can be made to fit the observed distribu- 
tion of HVC velocities, many of them also predict some fraction of clouds indistinguishable 
in velocity from clouds in the Galactic disk, with |li,sR| < 90 km s~^ (e.g. figure 3 in Wakker 



1991): we call these low- velocity halo clouds (LVHCs). In particular, recent numerical sim- 



ulations have shown that if HVCs are formed in a cooling Galactic halo scenario, we would 



expect to see an equal number of LVHCs formed in this process as HVCs ( Sommer-Larsen 



2006, Peek et al. 2008, hereafter PPS-L08). Indeed, the ratio of LVHCs to HVCs could be 



a powerful metric in distinguishing between these different HVC models. 



1.1. Models that Predict LVHCs 

As an example of the power of LVHCs to discriminate between models, we compared 
two very simplified models: one of clouds that inherit their velocity from the Galactic disk 
(the 'disk-associated' model) and one of clouds that inherit their velocity from a static halo 
(the 'halo-associated' model). These models in no way describe the physics involved in the 
production and evolution of LVHCs and HVCs, but rather are intended to highlight the 
effects of kinematically differing populations on LVHC and HVC statistics. 

The 'halo-associated' model is one in which clouds randomly populate a Galacto-centric 
sphere of radius 50 kpc, having a normal distribution of velocities parameterized by a velocity 
dispersion an and an overall infall velocity towards the Galactic Center of Vh- The 'disk- 
associated' model is one in which the clouds randomly populate a Galacto-centric cylinder 
with radius 15 kpc and \zmax\ = 10 kpc, with velocities inherited from the rotation of the disk 
{V = 2200 km s~^) nearest them, plus a random normal velocity distribution parameterized 
by (Jd, and an overall infall towards the disk of Vd- 

We then ran a Monte Carlo simulation in which we construct a random sample of clouds 
drawn from the distributions in each these models. We then 'observe' the distributions from 
the solar position {Rq = 8.5 kpc) to determine the histogram of HVCs as a function of Vlsr. 
We compare these histograms of HVCs to the histogram of HVCs in the updated Wakker and 



Van Woerden catalog (Wakker & van Woerden 1991) and thereby fit each of these models 



by allowing only those two parameters, V and a, plus an overall scaling for each model to 



vary. The results of this fitting process are shown in Table [T] 



Model name 


V [km s-i ] 


a [km s ^ ] 


^LVHc/^HVC 


Halo-associated 


-80 


79 


0.69 


Disk-associated 


11 


105 


1.58 


PPS-L08 simulation 


N/A 


N/A 


0.76 



Table 1: A table of the fitted parameters for each toy model and the corresponding ratio of 
LVHCs to HVCs. V = Vh.Vd for Halo and Disk, respectively, similarly, a = aH,crD- Also 
included is the LVHC-to-HVC ratio found in the simulations in PPS-L08. 

In Figure [T] we show the best fit of each of these two toy models. We also show the 
results of the fully cosmological simulation from PPS-L08, which is not fit to the data other 
than an overall scaling constant to take into account the limited resolution of the simulation. 
We see that these two toy models and one simulation each fit the observed distribution of 
HVCs reasonably well, but have more than a factor of two variation in the number of LVHCs. 
It is interesting to note the excess of observed clouds near the HVC cutoff velocity of 90 km 
s~^ , perhaps indicative of contamination by non-halo clouds. The variable LVHC result 



is consistent with plots in Wakker (1991) of toy models of HVCs stemming from accreting 
clouds and 'fountain' clouds. The broad, flat distribution in the 'halo-associated' model 
naturally stems from the 'smearing out' in the velocity domain by the solar motion in the 
GSR frame, and therefore leads to a much lower number of LVHCs. This wide variation in 
LVHC number across models, ranging from 69% to 158% of the HVCs observed, demonstrates 
the possibility that LHVCs can help to distinguish among HVC origins. 



1.2. The dust-to-gas ratio method 

One of the most important results from the Infrared Astronomical Satellite (IRAS) 
was the discovery of the strong correlation between the Galactic Hi column density and 



the infrared diffuse emission for |6| > 10° (Boulanger & Perault 1988). The conclusion 



reached by these authors, and others, is that there is a relatively uniform distribution of 
interstellar radiation field (ISRF) in the Galaxy and a relatively constant dust fraction in 
the neutral, diffuse Hi. The ISRF heats the dust, which reradiates in the infrared, and 
thus infrared flux is correlated with Hi column density. As Hi is typically optically thin at 
high Galactic latitudes, the column density is simply proportional to the observed velocity- 
integrated specific intenisty. The observed ratio of infrared emission to Hi column density is 
called the 'FIR-to-Hl ratio' or, put more simply, the 'dust-to-gas ratio' (DGR), and can vary 
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by factors of a few across the high-latitude sky (Boulanger & Perault 1988). (Note that in 
this work we use DGR to mean the ratio of the observed intensity of dust emission to that 
of HI gas, rather than the ratio of their masses, as it is sometimes used in the hterature.) 
It was also shown that HVCs do not follow this correlation; they have no detectable dust 



emission (Wakker & Boulanger 1986, hereafter WB86). WB86 speculated that this paucity 
of emission is because either the clouds do not have very much dust in them to start with, 
or because they are too far away to be bathed in the ISRF of the Galaxy. In either case, this 
lack of infrared flux is a velocity-independent observable that is associated with HVCs and 
therefore also with LVHCs. We can use this effect as a discriminant between low-velocity 
clouds that are physically part of the Galactic disk and LVHCs. 

The aim of this paper is to use the lack of infrared flux associated with HVCs, but not 
other clouds, to find the first LVHCs. We also place new limits on the DGR of known HVC 
complexes, present a new method to determine a more accurate DGR and determine the 
variability in the DGR of other lower velocity clouds and one dwarf galaxy. The paper is 
organized as follows: in ||2]we explain the observations in both infrared and Hi, in ^the data 
reduction methods for each are explained, in ^ we explain our technique for determining 
the contribution to the infrared flux from clouds, in ^we describe how we choose the clouds 
we investigate, we give our results in ^ and we conclude in ^ 



2. Observations 

The observations for this project were both done in the capacity of surveys: The IRAS 
survey provides maps of the Galactic, diffuse mid-infrared dust emission at high latitudes 
and the Galactic Arecibo L-band Feed Array survey in Hi (GALFA-Hl) provides maps of 
the Galactic Hi in the 21-cm hyperfine transition of Hydrogen in the same region. 



2.1. GALFA-Hi 

The radio data were obtained from the GALFA-Hl survey. ALFA is a 7-beam array 
of receivers mounted at the focal plane of the Arecibo 305m telescope Gregorian dome. 
GALFA-Hl uses this instrument to study Galactic Hi in the hyperfine transition of neutral 
hydrogen at 1420.405 MHz. The spectrometer, GALSPECT, has a velocity resolution of 
0.18 km s~^ (872 Hz) and a bandwidth of 1531 km/s (7.1 MHz) centered at km s~^ Local 
Standard of Rest (LSR) for each of the 7 beams and 2 polarizations. Each beam in the 
focal plane array has a FWHM of ~ 3.5', though the 6 non-central beams have significant 
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Fig. 1. — In solid black is a histogram of HVCs by LSR velocity in the updated Wakker and 
Van Woerden catalog (Wakker & van Woerden [1991 ) that fulfill the typical HVC criteria 
|^ls_r| > 90 km s~^ and \Vdev\ > 60 km s~^ (see Wakker 2004 for details). The vertical 
dashed lines show the region in which clouds are no longer considered HVCs by their LSR 
velocities. The best model fits (see ^1.1 and Table [ijfor details) are shown for the halo model 
(red, dotted) and for the disk model (green, long dashed). We also display the cosmological 
simulation data from PPS-L08 (blue, dot-dashed). 
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asymmetrical sidelobes towards the outside of the array sky-footprint displaced ~ 5'. The 
GALFA-Hl survey is conducted in both simple drift and basketweave (or 'meridian nodding') 
modes. Drift observations, taken commensally with the Arecibo Legacy Fast ALFA survey 
(Giovanelli et al. 2005), have a typical per-beam integration time of tint ~ 14.4 seconds per 
beam, whereas basketweave observations have tint ~ 4.8 seconds per beam. Some regions 
were observed with multiple passes in either mode, thus increasing the sensitivity of the 
observations. The final sky coverage of the survey will be ~ 13000 deg^, covering the entire 
right ascension range from —1° declination to 39° declination. At the time of writing less 
than 5000 deg^ of this region had been surveyed, collected between Spring 2005 and Spring 
2007 in the ongoing GALFA-Hl survey. 



2.2. IRAS 

IRAS observed the Galaxy over a period of 300 days in 1983 in the infrared wavelength 
bands centered on 12, 25, 60 and 100 fim ( |Beichman|p87| ). 98% of the sky was surveyed, al- 



lowing relatively complete overlap with the aforementioned 21-cm line observations. Though 
more sensitive and higher resolution observations have since been conducted with the Spitzer 
Space Telescope, the publicly available IRAS data were chosen as a comparison data set for 
their relative completeness and comparable resolution to GALFA-Hl. In this work we only 
concern ourselves with the 60 and 100 /im bands, which trace the cooler dust associated with 
Hi. 



3. Data Reduction 



3.1. GALFA-Hi 



To generate the Hi data cubes, raw data were reduced in the GALFA-Hl standard 



reduction pipeline (version 2.3), the details of which are described in Peek & Heiles (2008). 



See also Stanimirovic et al. (2006) for a description of the reduction process. Corrections 



are made for the IF bandpass, gain variation in the the receivers, impedance mismatches 
in the signal chain, static RF fixed-patten noise ('baseline ripple') and overall system gain. 
We reduce the contamination from all of these effects significantly, but the resulting maps 



are not wholly without systematic effects. In §5.3.1| we address the effects these systematics 
may have on our results. 



In any region of the sky under consideration, we define the 'cloud' and the 'zero-velocity' 
components of the gas. The cloud must be separated in velocity from the bulk of the 
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Hl gas, so we can make a simple cut between the two. This means that we cannot typically 
inspect clouds amongst the zero- velocity gas, which limits our sample considerably. We then 
integrate each component over their respective velocity ranges to make 2D maps of each (see 
Figure |2]). Any data containing 'glitches' are masked out. Low-amplitude striation can occur 
in the maps. In the case of the clouds it is typically dominated by residual baseline ripple 
effects not reduced by our initial data reduction. This ripple is effectively an additional 
spurious signal that varies slowly as a function of frequency. The striation is primarily along 
lines of constant declination, consistent with the drift pattern in which most of the data 
were taken. In this case, we reduce this effect by removing an average striation, measured 
in adjacent regions with the same declination. These adjacent regions contain no Hi flux, 
and so are good measures of the striation to be removed. The striations in the zero-velocity 
gas, also typically consistent with the drift pattern (constant declination), primarily stem 
from errors in our fits to the varying beam gains. The zero-velocity maps are cleaned by 
first averaging the map across declination. This declination-averaged map is then smoothed 
to a scale slightly larger than the striations. The smoothed component is then subtracted 
from the declination-averaged map to determine the striation component. This striation 
component is then in turn subtracted from the original zero-velocity map to yield a map 
with significantly fewer striations. 

We also note that the Arecibo telescope is known to have non-trivial distant sidelobes 
(stray radiation), wherein some radiation from other parts of the Galactic Hi sky can contam- 
inate the observations. Comparisons to the stray-radiation corrected Leiden- Argentina-Bonn 



survey (Kalberla et al. 2005) do not show detectable differences in the shape of the HI spec- 
tra, though we do allow that some stray radiation will certainly contaminate our data and 
lead to increased scatter in our results. 



3.2. IRAS 



The IRAS extended emission data were originally reduced in 1984 and 1986 and were 
released as the SkyFlux atlas. These data were later reprocessed with more knowledge of 
the IRAS instrumental systematics and released as the IRAS Sky Survey Atlas (ISSA) in 
1991 and 1992. This reduction, though exceedingly useful, was of limited value because it 
suffered from significant striping, 'glitches,' zero level effects and zodiacal light. The data 



were once again re-reduced by Schlegel et al. (1998) to build an all-sky extinction map. 
This significantly reduced the striping and 'glitching' in the data. Finally, the data were 
re-reduced again by Miville-Deschenes & Lagache ( 2006[ ) (hereafter M-DL06), increasing the 
effective resolution in both the 60 /xm and 100 /xm bands to ~ 4' (comparable to that of the 
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Fig. 2. — These plots show a typical region of the sky we examine and the associated cloud 
and zero-velocity gas column density maps. The top plot is the average spectrum over lines 
of sight where we detect the L8 cloud in the Hi data. The light gray velocity range is the 
range over which we integrate to determine the zero-velocity gas column density; the dark 
gray is the range for the cloud. Note that we include all Galactic Hi along the line of sight not 
included in the cloud to determine the zero-velocity gas column, even if this includes some 
clouds that are separable from the main Galactic emission. The two maps are of the cloud 
column density (middle) and the zero- velocity column density (bottom). The checkerboard 
area represents regions that are either unobserved by GALFA-Hl or the data are omitted 
due to issues of data quality. 
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GALFA-Hl survey) as well as fixing significant issues with zero-point offsets and zodiacal 

light. Wc Tisc this reduction, called IRIS, because of the improved angular fidelity in the 
bands of interest. In addition to the reduction outlined in M-DL06, we search the data for 
point sources down to a limiting magnitude of 1 Jy, as we are not interested in the infrared 
emission from stars. We remove these point sources by fitting with a typical PSF as quoted 
in M-DL06. The PSF of IRAS is not simply Gaussian, so the brightest stars (typically < 1 
/ field) are difficult to remove and are therefore masked out of the final comparison 60 /im 
and 100 /im images. 

4. Tcirget Selection 

We are interested in constraining the DGR for several different types of objects: HVCs, 
lower-velocity clouds, and Local Group dwarf galaxies. Note that the range of clouds with 
I^lsrI < 90 km s"-*^ has historically been divided into clouds of low velocity and clouds of 
intermediate velocity. As this distinction is irrelevant to this work, we will hereafter refer 
to all clouds distinguishable from the main Galactic zero- velocity emission with |Vlsr| < 90 
km s"-*^ as low-velocity clouds, or LVCs. The DGR of HVCs can be used to constrain their 
origin and confirm previous studies. The DGR of LVCs can be used to discern their possible 
membership the the Galactic halo and similarity to HVCs. The DGR of Local Group dwarfs 
may allow us to differentiate them from Galactic gas. Our samples of these objects are 
described below. 



4.1. HVCs 

Our first goal is to inspect any HVC complexes with members observed in the current 
GALFA-Hl data set and measure their dust-to-gas ratio. It is important to do this for two 
reasons. First, these new data sets and techniques yield higher fidelity assessments of the 
HVC DGR than have been previously made for the available HVCs, which may place tighter 
limits on the physical conditions of HVCs. It is also possible that we may detect dust in 
these HVCs. Note that a detection of dust in HVCs, if it is at a very low level as expected, 
will not derail our fine of reasoning for using the DGR to determine whether specific LVCs 
are actually HVC analogs: as long all HVCs are distinguishable from normal Galactic gas 
via their DGR, the reasoning holds. Secondly, it is important to establish reasonable upper 
limits for the DGR of HVCs to compare to our LVC sample; indeed, the main thrust of this 
work is to find LVCs that are similar to HVCs in their measured DGR, so a baseline upper 
limit for HVCs is crucial. 
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To find HVCs in our map we cross reference all the clouds in the WVW91 catalog with 
the coverage map of the current GALFA-Hl data set. We then find the highest fiux cloud from 
each of the HVC complexes in the GALFA-Hl data set and perform image cleaning methods 
described in [3A on that cloud's region of the sky. In many cases we find that individual 
clouds in WVW91 can be resolved into multiple components, an unsurprising result given 
the superior resolution of the GALFA-Hl data set. For a given HVC complex, the region we 
use typically includes the bulk of the fiux observed within the GALFA-Hl survey to date, 
but not necessarily the majority of the total fiux of the complex. For instance, only the 'tail' 
of complex C is visible to Arecibo, so our measurements are only of this region of complex 
C. We apply the 'displacement map' method described in §5.2| to the resulting maps to 
determined values of DGR for each cloud and associated errors. 



4.2. LVCs 

The method for finding LVCs is less straightforward than for finding HVCs, as we do 
not have a catalog. While some LVCs are cataloged in HVC catalogs (e.g. WVW91), it 
would be impossible to determine from such a lower-resolution Hi data set which are clouds 
separate from the zero- velocity or 'disk' gas and which are intermingled with the disk gas. 
Given this situation, we search through the entirety of the relatively contiguous portions of 
the GALFA-Hl data set (~ 4000°) by eye to determine likely targets for our method. These 
targets must have |Vlsr| < 90 km s~^ but be distinct from the disk. We allow some very 
small portion of the fiux to overlap from cloud to zero-velocity gas, as long it is in the far 
wing of the clouds, and thus does not significantly contaminate our column density maps. 

We found a total of 9 clouds for which dust measurements could be made, ranging 
from \b\ = 84 to \b\ = 13 and with fluxes ranging from 170 to 6600 Jy km s~^ (see Table 
[2]). A similar number of clouds were rejected for various reasons. These reasons include an 
unconstrained value of -Da,C) stemming from a paucity of flux or extreme variability in the 
zero-velocity DGR; signiflcant systematic errors in the GALFA-Hl data set or contamination 
in the IRIS data set; and mixing between the cloud and zero- velocity gas, thwarting the 
construction of accurate column density maps. The typical LVCs we flnd are relatively 
compact, clumpy clouds like HVCs, and are quite visible when examining the Hi data cubes. 
We note that this sample of clouds is not necessarily statistically representative of all LVCs. 
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4.3. The Leo T Dwarf 



In addition to targeting Galactic clouds, the main thrust of this paper, we also investi- 
gated using this method to constrain the DGR of Local Group dwarf galaxies and potentially 
find new ones by examining the newly-discovered gas-rich dwarf galaxy Leo T. Gas-rich dwarf 
galaxies around the Milky Way are very difficult to distinguish in Hi from compact Galactic 
clouds. Leo T, for example, has a radial velocity of only 30 km s^^ with maximum HI in- 
tensity of 3.5K, comparable to typical Galactic local gas clouds (although relatively distinct 
from HVCs or LVHCs, which have a maximum HI intensity of ~ 1 K, and are typically less 
compact). Finding such galaxies via their HI signature alone is of great interest, as recent 
studies have shown that dwarf galaxies beyond 300 kpc can have very high gas fractions. 



and therefore may be undetectable in starlight (Grevich & Putman 2008). Observations 



compiled by Lisenfeld & Ferrara ( 1998 ) show that dwarf irregular galaxies with HI masses 
below 1O^-^M0 have very low or undetectable DGRs. Thus, if Leo T, the only such confirmed 
dwarf galaxy in our observed area, does indeed have very low DGR, it is a confirmation that 
the same method by which we are attempting to distinguish LVHCs from other LVCs may 
be used to distinguish other dwarf galaxy candidates from compact Galactic HI clouds. 



5. Methods 

Since the Hi column densities of HVCs (and presumably the analogous LVHCs) are 
typically no more than 10% of the high-latitude Galactic Hi column density, we expect that 
even in the case of a high dust-to-gas ratio cloud the infrared contribution from the cloud 
will be quite limited. Therefore, the methods by which we determine the dust-to-gas ratios 
of these clouds require significant attention to avoid contamination from systematics. 



5.1. The simple approach 

To determine the DGR of a cloud, the simplest method is to isolate the column densities 
from the cloud over an area of sky, Nq (a, 5), and from the zero- velocity gas, Nz (a, 5), (see 



3.1) and solve the set of linear equations 

h («, 5) = D^^cNc («, 5) + D^^zNz (a, 6) + if^, (1) 

to minimize the sum of the squares of the errors, where I\ (a, 5) is the infrared intensity at 
wavelength A over the region and K is an arbitrary offset in the infrared data. A least-squares 
solution for the equation for the range of values of a, 5 across the maps is implemented to 



-VI- 



solve for -Da,c and D\ z-, the values of the DGRs of the cloud and zero- velocity components 
respectively. Note that l\ (a, S) can stem from other sources, such as dust associated with 



molecular gas (e. g. Gillmon & ShuU 2006, Gillmon et al. 2006) and dust associated with 



ionized gas (although see Odegard et al. 2007). These contributions are not considered in 



this work, and are therefore sources of error to our fits. 

We denote the 'true' values of the DGRs as D\fj and -D^^ and the 'true' offset as 
in contrast to the above measured values. This technique is dependent upon the assumption 
that D\ 2, a^iid K\ each have a single value in the region independent of position. If the 
zero- velocity DGR were to have significant spatial dependence within the region (D\ z ~ 
D\ 2, (c^! ^T)-, Dx,c would depend critically on whether regions of high Nc {a, 6) overlapped 
regions of relatively low or high D*^ ^ (a, 5) as compared to the measured -Da,z- Again, it is 
important the realize that in this method variations in the DGR can include both variations 
in the true dust emissivity or fraction and variations in the contribution from other sources 
(molecules, ionized gas) to I\. 

This effect is illustrated in Figures [3} |4| and [5j In Figure |3] we show the integrated 
column density and infrared intensity from a small region of high latitude sky, column 
density Nz {a, 6) and infrared intensity Jioo («, 5) respectively. We also show the integrated 
column density of a typical cloud from another region of the sky, Nc {a, 5). Figure |4] is a plot 
of Nz vs Jioo, with a line fit to the data. It is clear from this plot that there is significant 
scatter from a simple linear fit. This scatter has RMS noise ~ 0.15 MJy sr~^, meaning that 
for a cloud of typical column density lO^^cm"^ to have infrared flux equal to the scatter, it 
would have to have -Da,c — 150 x 10~^^MJy/sr cm^. This residual scatter is mapped in the 
left panel of Figure [5} demonstrating our claim that there is significant spatially correlated 
variability in ^. The right panel of Figure [5] is a demonstration of the variability induced 
in Dx c by moving the cloud shown in the right panel of Figure |3] around in a and S. There is 
no infrared flux associated with the cloud, but as the cloud is positioned in different locations 
in the sky, the -Dioo.c solution to Equation [T] varies dramatically. The right panel of Figure |5] 



is an image of this variation (see ^5.2 and Equation [2] for details), which is similar in form to 
the left panel of Figure [s] Indeed, -Dioo.c varies from —75 to 130 x 10~^^MJy/sr cm^, while 
Dioo^z — 35 X 10~^^MJy/sr cm^. This shows that the noise in such a measurement typically 
exceeds the value we are trying to measure. 

We draw the somewhat unfortunate conclusion that the application of Equation [T] cannot 
yield a believable -Da,c result in the case of the majority of Galactic clouds. This conclusion 
is not terribly surprising, as we know that Equation [T] does not include other sources of 
FIR radiation, and so can have significant, spatially coherent contaminants. If the cloud 
happens to be along the line-of-sight towards a region with particularly consistent Dx,z and 
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the cloud has a very high relative column density [Nc/Nz ~ 1), it may be possible to use this 
method effectively, but for the majority of low column density clouds a new method must 
be developed. In the case of the claim of detection of dust in complex C using an equation 
equivalent to eqn. [l]by Miville-Deschenes et al. (2005), in contrast, we find complex C does 
not have a very high relative column density. Our analysis therefore casts significant doubt 
on this detection. 



5.2. The displacement map method 

In the previous section we showed that displacing a cloud relative to the Galactic emis- 
sion can change the measured relationship between the Hi and IR flux, the measured dust- 
to-gas ratio. We use this idea of displacing clouds relative to the Galactic and IR emission 
to develop our new measurement technique. We call this the 'displacement map' method. 

The displacement-map method is a member of a broad class of methods called 'tem- 
plate matching', developed by the computer vision community over the last few decades 



(e.g. Rosenfeld 1969). In a template matching method, the goal is to determine whether, 
where, and to what extent an image, /, is represented inside another image, g, or to deter- 
mine if / and g are similar. This is done by using a 'template' of / and applying a wide 
variety of algorithms to determine the significance of a match to g. In computer vision 
the position, scale, and rotation of the template of / (among other possible modifiers) are 
generally unknown, thus making this problem complex and computationally intensive. 

In the displacement-map method we confront a similar problem - we expect to see an IR 
contribution from a cloud in our infrared map, but the map contains other signals and noise 
as well. We expect that the IR contribution to the map from the cloud should look like the 
cloud (i.e. we assume Dc is a constant); this is the template we use in the displacement-map 
method. As opposed to most template matching problems in the field of computer vision, 
we are only interested in how well the template matches at the position of the cloud on the 
sky. What follows is a qualitative description of the method. 

We quantify how well our template, the cloud, matches each position in the residual IR 
map after we have fitted out the dominant contribution from the zero-velocity gas. This is 
our displacement map, Dc (Aa, A6) (note that we drop the subscript A in this section as 
the method applies to either FIR band and to reduce clutter). At the center of this map 
will be our 'signal', contaminated by chance correlation between the residual scatter and the 
template. To reduce the effect of these fluctuations, we filter out large-scale structures from 
the displacement map - this highlights the relatively sharp peak we expect at the center of 
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Fig. 3. — The left image is the disk Hi column density from a region at high Galactic latitude 
from GALFA-Hl (1 = 200, b = 50), while the center image is the 100 /im emission in the same 
region of sky from the IRIS data. We specifically choose a region with no gas beyond the 
zero-velocity disk gas, so there is no cloud-associated infrared-flux i.e. -D^^; = 0. The right 
image is the integrated Hi column density from GALFA-Hl of a typical cloud we investigate 
from a different region of sky. We use these three maps as examples of Nz (a, S), /loo («, 
and Nc{a,5), respectively to demonstrate the difficulty in using Equation [l] to determine 
Dx,c- 
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Fig. 4. — This is a plot of the Nz (a, 5) vs. /loo («, 5) data from Figure|3| with a superimposed 
hnear fit to the data. 
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Fig. 5. — The left image is of the residuals of the line fit shown in Figure |4| Jioo (a, 5) — 
KioQ — DiQQ^zNz {a,6). It is clear that there is significant spatial correlation of residuals 
across the region - this is typical over the regions we have examined. The right plot is the 
Dx,c (Aa, A5) 'displacement map' (see ^5.2 and Equation [2| for details). It shows the effect 



of placing the cloud at different locations throughout the image (Aa, AS) on the calculated 
Dx,c- Therefore, the variability in this map shows the dependence of the measurement of 
Dx c on the region along the line of sight to the cloud. The visual similarity between the 
left panel and the right panel demonstrates the effect that relatively 'over-dusty' regions in 
the zero-velocity gas have the strong tendency to drive up the Dx^c- 
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the map. We then measure the amphtude at the center of this filtered map and the noise in 
that map. These measurements yield a gauge of the dust emission from the cloud and the 
accuracy of that result, respectively. 

Here we present the detailed, rigorous description of the method. We solve the equation 

I (a, 6) = Dc (Aa, A6) Nc {a -Aa,6- A6) + Dz (A«, A6) Nz (a, 6) + K {Aa, A6) , (2) 

for a grid of displacement values of Aa and A6, developing the displacement maps K (Aa, A6), 
Dc {Aa, A6) and Dz {Aa, A6). We expect that there should be little variation across 
K {Aa, A6) and Dz {Aa, A6), as the values of these fits are primarily dependent upon the 
correlation between Nz {a, 6) and I {a, 6), which are independent of the variation of position 
of the displaced cloud on the sky, Nc {a — Aa, 6 — Aa). By contrast, Dc {Aa, A6) should 
have strong variations depending upon the strength of any cloud infrared emission signal in 
/ {a, 6) and the uncompensated spatial variations in the values of K* and D*^ (as shown in 
Figure |5]) . 

Our new displacement map technique depends upon the fact that if a cloud does indeed 
have associated dust flux it will have a predictable effect upon the Dc {Aa, A6) displacement 
map. We take the infrared emission to be decomposed into two parts: 

I{a,6)=I^{a,6) + I^{a,6) (3) 

where 

l''{a,6)=D*cNc{a,6) (4) 

and {a, 6) is the component of the flux associated with the low velocity gas and residuals. 
Note that we use the superscript C or Z to denote that the quantity is associated with 
infrared flux from either the cloud or zero- velocity component, whereas the subscript C or Z 
is used to denote whether the quantity is associated with the Hi column from either the cloud 
or zero- velocity component. As an example, Dc {Aa, A6) represents the entire displacement 
map of interest, while Dc {Aa, A6) is the component of the displacement map associated 
with the flux that comes from the cloud itself, This construction allows us to determine 
this influence of cloud-associated dust emission on Dc {Aa, A6). 

Since we solve Equation |2] with a linear least-squares fit, the solution can be decomposed 
into two solutions, one for each of the components of / (a, 6). We are interested in the solution 
to the J*^ component of the equation: 

I^' {a,6) - {Aa,A6) = D^ {Aa, A6) Nc {a - Aa,6 - A6) (5) 

+Dg {Aa,A5)Nz{a,6). 
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Since to first approximation Nz (a, 5) should be uncorrelated to (a, 5) in general (D^ = 0) 
and (a, 5) has no 'offset' component, Equation [s] can be reduced to 

(a, 5) ~ (Aa, A5) A^c (" - Aa, 5 - A5) (6) 

Substituting in Equation |4] we find 

D^A/'c (a, 5) = Dg (Aa, A5) A^c (a - Aa, 5 - A5) . (7) 

This is equivalent to a simple least-squares fit of slope with y- intercept = 0, whose solution 
is simply 

Dg (Aa, A(5) = Dl 



.c,A. A x^ _ E A^c ^^c {a - Aa,6 - A6) 



The ratio of the sums is simply the normalized two-dimensional auto-correlation function of 

In practice, i?„ {Nc (a, S)) varies much more rapidly near the origin than the overall variabil- 
ity of Dc (Aa, A6), which typically varies relatively smoothly (e.g. the right plot in Figure 
[5|. We can remove this smooth variability from Dc (Aa, A5) by subtracting a smoothed 
(i.e. convolved) version of the displacement map: 

U {Dc (Aa, A6)) = Dc (Aa, A6) - k (Aa, A6) * Dc (Aa, A6) , (10) 

where k {Aa, A6) is our Gaussian smoothing kernel. This is equivalent to either applying an 
'unsharp mask' or filtering out low spatial frequency modes in the Fourier domain. While 
this unsharp masking will have the effect of diminishing the contribution from spurious large- 
scale fluctuations in the displacement map, it also somewhat diminishes the strength of the 
signal of interest, Z)g(Aa,A5), at the origin as well. Since unsharp masking is a linear 
operation, 

U {Dc {Aa, A6)) = U (Dg (Aa, A6)) + U (Dg (Aa, A5)) . (11) 
Using Equation |8] we then solve for D^: 

„ . . ^-U(D^ (Aa. AS)) L „ , . „ 

(12) 



U {Dc {Aa, AS)) U„=o,a5=o - U {D^ {Aa, A5)) U,=o,a5=o 



U{Rn{Nc{a,6))) 



Aa=0,A(5=0 



All of the right side of this equation is known except U {^Dq {Aa, AS)), which we ignore as 
small and take to be contribution a to the error. 
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To determine the error in this measurement, we simply calculate the standard deviation 
(RMS noise) in this scaled, unsharp-masked displacement map away from the origin: 



^ U{DciAa,AS)) ^ 



U{Rn {Nc{a,5))) 



Aa=0,A<5=0 



(13) 



We set the size of the Gaussian smoothing kernel, k (Aa, A5), to a FWHM = 10'. This 
size is chosen simply to maximize the signal-to-noise. This displacement map method allows 
us to take into account accurately the effect of the errors in our fits to the zero-velocity gas 
without attempting to fully parameterize them. 

The displacement map method is illustrated in Figure [6] The six plots in Figure |6] are 
all on the same scale of one arcminute per pixel. The upper left plot is a reproduction of the 
right plot from Figure [sj Dc (Aa, A5) in Equation |2| shown here for direct comparison to the 
next plot. The plot in the upper right is similar to the previous plot but of Dq (Aa, A(5) as in 
Equationpl with D*(j 7^ by construction. In this case we set = 36 x 10^^^ MJy sr~^ cm^. 



The maps are very similar, except a very slight increase in intensity at the center of the 
map, corresponding to the detection of the cloud dust. The middle left plot is the difference 
between the two upper plots, which is to say Dq (Aa, A(5) — (Aa, A(5) = (Aa, A(5). 
The middle right plot is of Rn{Nc {a-i^)). It is empirically evident that the assumptions 
required to state Equation [6] are justified in this D'^ (Aa, A(5) is extremely similar 

to Rn {Nc {a, S)). The bottom left plot is of U {Dc (Aa, A6)). This shows how the central 
spike that comes from (Aa, A6) can be measured in the noisy background. The right 
plot is of U {Rn {Nc {a, 6))), which we compare to the bottom left plot to determine Dc- 



5.3. Tests of the displacement map method 

Here we test the efficacy of the displacement-map method and the scale of a few possible 
systematic contaminants. 



5.3.1. Measurement accuracy 



To test the claim in [5.2 that we can accurately measure the DGR in clouds (as well as 



the error in that value), we construct a numerical experiment. As in Figures [3] and [6| we find 
an arbitrary high-latitude region of the sky containing no clouds and measure the Nz (a, 6), as 
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Fig. 6. — Six plots on the same scale of one arcminute per pixel, demonstrating the displace- 
ment map method. See the end of ^5.2 for details. 
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well as [a, 6). We then find a cloud from a different region of sky and use it as our Nq [a, S) 
and assume a (j. We construct the dust map, Ix (a, 5) = (a, 6) + D\ ^Nc {a, 6). We 
then follow the displacement map method to measure -Da,c and the associated errors. We 
repeated this experiment with three different regions of the sky, three different measured 
clouds, each in 16 different positions, using 4 different values for D^qq,^: 15, 45, 75, and 150 
xlO~^^ MJy sr~^ cm^. We find using these 384 experiments that the values predicted by the 
method are consistent with the input value of D*^ ^ to better than 15% of the predicted 1-cr 
errors and that the predicted l-cr errors themselves are accurate to better than 5%. 



5.3.2. Gain and baseline ripple effects 



In section we described our methods for reducing the effects of baseline ripple and 
gain fluctuations. While these methods are successful in dramatically reducing these instru- 
mental effects, the residuals of both of these contaminants can produce visible effects in the 
zero-velocity and cloud maps. In both the case of residual gain fluctuations and residual 
baseline ripple, small correlations are induced between the measured column density of the 
zero-velocity map and the cloud map. This correlation can have the consequence of reducing 
the measured Dx^c- 

We test these effects by conducting numerical experiments similar to those conducted 



in { 5.3.1 We find a baseline result by using the displacement map method to determine 
Dx^c for both the Ix («, 5) and the Jf (a, 5) dust map. The results are consistent with the 
inputs, i.e. the derived -Da,c for the Jf (a, 5) map is consistent with zero and the derived 
Dxfi for the Ix (a, 5) is consistent with D\ (j. While the Nz (a, 5) and Nc (a, 5) are already 
contaminated by these effects, they are uncorrelated, and so should only add to the noise, 
rather than biasing the results of the fit. To add in this correlated noise, we then examine 
another region of sky for typical baseline ripple effects and typical gain effects. We apply each 
of these effects independently to the Hi column density maps, and re-run the displacement 
map technique for each set of contaminations. In each case we see an effect on the final 
Dx^c of — 5%, in the expected sense of decreasing the effect of the fit. Since it is difficult to 
determine the exact amplitude of this effect in the case of each cloud we consider, we add a 
10% systematic error to account for this effect. 
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5.3.3. Fit region effects 

In previous sections we have referred to the 'region' over which the equations are eval- 
uated, spanning some range in a and 5. We now ask how sensitive our solution for -Da,c 
is to variation in the size of that region. An experiment similar to the above section was 
conducted, wherein instead of adding noise to the maps, we vary the area of the region 
under consideration, keeping the size of the cloud fixed. We find that varying the area of 
the regions by factors of 3 has no detectable effect on the accuracy of the method. 



5.3.4. Effects of fitting routines 

The standard least-squares or chi-square fit only takes into account errors in one variable 
in the problem. To take fully into account errors in all dimensions, one can employ Jefferys' 



method, an iterative chi-square technique (Jefferys 1980). Unfortunately, this method can 



be somewhat unstable and is strongly sensitive to the assumed errors and covariances, which 
may be difficult to determine robustly. The method can also be very computationally inten- 
sive, particularly when being employed many times as described in §5.2 For these reasons we 



wish to determine the sensitivity of our results to the methods by which we solve our linear 
equations, and specifically whether we get accurate results using a standard least-square fit. 
We performed tests similar to those in the previous sections, but varied our fitting method, 
alternately using Ix (a, 5) or Nz (a, 5) as our independent variable, assumed to have no errors. 



While the values measured by the 'simple approach' (^5.1) did vary significantly between 
fixing different independent variables, and did not match our input D\(j, our results using 
the displacement map method were effectively identical to each other and matched D\ ^ to 
within the measured noise. Having empirically shown that the more complex multi-error 
methods are unnecessary, we use the standard least-squares fit method to determine Dx c- 



6. Results 
6.1. HVCs 

The results for HVCs, as shown in Figure [7] and Table |2| primarily confirm the result 
from WB86 — HVCs have no detectable dust emission. Some of these HVCs (WA, G, and 
P) have relatively small errors, clearly setting them apart from the typical measured DGR in 
the Galactic disk, Da.go ~ 20 x lO'^^MJy sr'^ cm^ and Da.ioo ~ 100 x lO'^^MJy sr"^ cm^ 
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LI 


84 
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20 


5.9 X 10^ 


39 ±20 


12 ±5.6 


0.81 


L2 
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60 


48 


6.6 X 10^ 


17 ±4.4 


6.4 ±2.0 


0.81 


L3 


322 


68 


30 


3.0 X 10^ 


172 ± 16 


64 ± 7.8 


0.86 


L4 


5 


53 


-60 


2.3 X 10^ 


153 ± 19 


45 ± 7.3 


0.86 


L5 


253 


60 


53 


3.1 X 10^ 


7.2 ±9.1 


5.9 ±4.6 


0.40 


L6 


32 


84 


-58 


1.7 X 10^ 


2.8 ±6.3 


1.9 ±3.2 


0.59 


L7 


47 


39 


-60 


2.7 X 10^ 


67± 16 


27 ±4.0 


0.91 


L8 
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-29 


43 


2.9 X 10^ 


27± 18 


21 ± 7.1 


0.72 


L9 


38 


65 


-85 


2.4 X 10^ 


34 ± 4.0 


11 ±1.7 


0.91 


Dwarfs 
















Leo T 


211 


13 


30 


7 X 10" 


8.5 ± 13 


-6.8 ± 5.1 


0.54 



^AU dust-to-gas ratio values are listed in 10 ^^MJysr ^ cm^. 

Note. — A table of objects observed in our preliminary survey of HVCs, LVCs, and a 
dwarf Galaxy, p is the measured correlation coefficient between the errors in the 60 and 
100 //m DGRs. Errors are 1-a values. 
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Fig. 7. — Two plots of the results of our analysis for 6 HVCs and 9 LVCs and the Leo T 
dwarf. The top plot focuses on the LVCs, the bottom plot focuses on the HVCs and the 
Leo T dwarf. HVCs are labeled by the complex to which they belong. Our data are shown 
by 1-0" ellipses, which take into account the correlation, p, in errors between the 60/im and 
lOOyum DGR measurements. In the top plot we have overplotted data from BP88 - this is the 
DOR for the total Hi column from 3 regions of the sky (Auriga, Orion and Lupus) as well 
as the Northern and Southern Galactic Caps beyond |6| = 50, the sky near |6| = 30 and the 
fit to an assumed csc(b) distribution. We also show the expected zero infrared emission for 
HVCs with a cross in the bottom plot. Due to lack of available DIRBE calibration the BP88 
data may have systematic errors at the 20% level. Note the detection of the dust in complex 
M, the non-detection of dust in L5 and L6, as well as the increased DGReo/im/DGRioo/im in 
LVCs compared to the total column results. 



( Boulanger fc Perault|1988 ). The results are less clear for the HVCs we sampled in complexes 
C and AC, which have larger error bars. Note that our measurements do not strictly rule out 
the values found in Miville-Deschenes et al. (2005). Also note that a few of the HVCs show 
negative DGR at low significance (< 2(t), which we take to be measurement noise. The most 
surprising result is a > So" detection of dust in a cloud in complex M. It seems that complex M 
does indeed have some small amount of heated dust. Complex M is a relatively low- velocity 
HVC complex and is on the border of the HVC/IVC classification. Indeed, some think that 
the cloud may be related to the IVC arch at a distance of ~ 1 kpc (van Woerden & Wakker 



2004, Wakker 2001). Complex M has been shown to be within 4 kpc by UV absorption line 
measurements (Danly et al. 1993) and Ca II K absorption (e.g. Keenan et al. 1995), but 



does not have significant lower limits. If it were indeed related to the IVC arch, that might 
explain the detection of dust, given the strong expected UV flux near the plane of the Galaxy 
(e.g. Putman et al. 2003 Bland-Hawthorn & Maloney 1999) and known solar metallicity of 
the IVC Arch (Richter et al. 2001). There are no absorption- line metallicity measurements 
in the literature as there has not been a high-resolution measurement of Hi towards the stars 
detected behind complex M. Measurements of [S II] and Ha by Tufte et al. (1998) point 
to a metallicity near solar, but there are significant errors in this measurement. A solar 
metallicity would be consistent with a Galactic origin for complex M. All other measured 
HVCs (WA, G, P and AC) are consistent with low dust dust content and/or negligible ISRF 
at the cloud. 



6.2. LVCs 

Results for the nine LVCs measured are also shown on the top of Figure [T] (ellipses) and 
Table |2} Our measurements of LVCs span a very large range in both DGRioo and DGReo- 
We find 2 such clouds, L5 and L6, which are consistent with zero DGR and inconsistent 
with normal disk gas. These are candidate LVHCs. We also find one cloud (L2) with DGR 
consistent with the DGR found in the complex M detection. L2 may be consistent with 
whatever kind of object complex M is (a nearby HVC or a fast-moving disk cloud), or it 
may be on the far end of the typical range of LVC DGRs. The remainder of the LVCs 
are distributed over a wide range of DGR values. These values roughly span the same 



range shown in the analysis by Boulanger & Perault (1988) of the DGR of various areas 



of the high galactic latitude sky, including star-forming regions, the Galactic caps and the 
overall sky averages beyond b= |30|: 10 <DGR6o < 55, 50 <DGRioo < 240, in units of 
10-22MJy/sr cm^. 

Of note is that the measured DGRs in these LVCs have a typically higher DGRgo/DGRioo 
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value than the measurements from Boulanger & Perault (1988). There are a number of possi- 



ble reasons for this effect. The effect may arise from an overall smaller grain size distribution, 
which would result in a higher temperature for the grains and thus 'bluer' grains. There 
also may be a larger number of very-small grains (VSGs); VSGs are stochastically heated 
by UV photons and thereby contribute significantly to the GOfim grain emission. Given the 



conclusion in Fitzpatrick (1996) that LVCs have less material locked up in grains due to 



grain destruction from shocks, the decrease in average grain size (or increase in VSGs) may 
indeed be the correct interpretation of these results. This result is consistent with previous 



work on the subject by Heiles et al. (1988), but inconsistent with observations by Deul & 
Burton (1990), who find a lower DGRgo/DGRioo value in such clouds. 



We find that it is very difficult to make statistically sound claims about the total number 
of LVHCs that exist based on these data, primarily due to the complex selection effects 
involved. We might argue that if we take the number of LVHCs observed and divide out 
both our detection efficiency for LVHCs and the area of the sky surveyed we might be able 
to extrapolate to the total number of LVHCs in the sky. There are a few difficulties in 
pursuing this line of reasoning. The first is that we detect a very small number of LVHCs. 
With only two clear cases of candidate LVHCs, we have very large error bars simply from 
Poisson statistics. It is also difficult to determine how many 'clouds' each LVHC comprises; 
to compare to the the WVW91 catalog we need to have the same criteria for determining 
cloud identity. We also encounter the problem of not having a very accurate measurement 
of our detection efficiency - certainly LVHCs at Vlsr = are undetectable to this survey, 
but at what velocity they become detectable varies significantly over the area of the sky 
observed. The detection efficiency is also strongly dependent upon the quality and depth 
of the GALFA-Hl data and IRAS data, which vary significantly. It is also important to 
note that since the area of the sky surveyed is rather small and that HVCs are known to be 
significantly clustered on the sky, any extrapolation may be strongly influenced by whether 
the area of the sky covered happens to intersect an LVHC complex. We therefore leave a 
statistical analysis of the LVHC population to a future paper with a larger cloud sample and 
quantiflable selection criteria. 



6.3. Leo T 



We flnd that Leo T has no detectable dust emission in either 60yum or 100/im, consistent 
with our conjecture that a gas-rich dwarf galaxy near the Milky Way will have cooler, or 
less, dust than the Milky Way itself. This flnding is also consistent with the analysis in 



Lisenfeld & Ferrara (1998), who show that gas-rich dwarf galaxies have signiflcantly lower 
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dust fractions (by mass) that the Milky Way. This result leads us to believe that the DGR 
method may also be useful for determining whether a compact HI cloud is simply a cloud in 
the Galactic disk, or whether it is a gas-rich dwarf Galaxy candidate. 



7. Conclusion 

We have shown that the GALFA-Hl data and IRIS reduction of the IRAS data make 
an excellent match for use in determining the DGR of neutral clouds in our Galaxy {^^■ 
We showed that the simple least-squares (or chi-squared) fit method for determining the 
dust emission from neutral clouds with a distinct velocity (LVCs & HVCs) is flawed in 
circumstances where FIR fluxes are low compared to that of the zero- velocity gas ( ^5.1 ). We 



developed the displacement-map method for determining the DGR of such neutral clouds and 
showed it to be an accurate and robust technique ( ^5.2 ^5.3). We applied this method to the 



9 LVCs and 6 HVCs found in the GALFA-Hl survey area mapped to date along with the Leo 
T dwarf galaxy (Q. Only one of the 6 HVCs, a cloud from complex M, showed a significant 
detection of dust. This cast into doubt the distance to, and history of, complex M, which 
may be related to more nearby, disk-associated clouds such as the Intermediate- Velocity 



Arch (^6.1). Two lower- velocity clouds, L5 and L6, were shown to have no detectable dust 
content, and we therefore consider them candidate LVHCs ( ^6.2[ ). We also find that most 
LVCs and IVCs have significantly hotter dust than typical Milky Way, as measured by their 
DGReo/DGRioo value, which we find to be consistent with observations of decreased dust 
depletion in LVCs. We also find that Leo T has no detectable dust emission, as we expect 
for a dwarf galaxy with relatively low ISRF (^6.3l). 



The clearest follow-up work needed here is to examine larger swathes of the Arecibo- 
accessible sky for HVCs and LVCs. As the GALFA-Hl survey accumulates data in the coming 
years we expect to measure many more clouds using this method to determine if there are 
more LVHC candidates, if any other HVCs show clear signs of dust, and if the trend of hotter 
LVC dust remains. If a large enough sample of LVC measurements can be built, we may 
be able to statistically address the question of how many LVHCs there are as compared to 
HVCs, thus giving us a new discriminant between HVC models (see Figure [T]). In particular, 
we are pursuing automated methods of searching for LVCs in the GALFA-Hl survey, which 
will allow us to determine our selection biases much more accurately and thus extrapolate 
to the full LVHC population. These cloud isolation methods are already showing significant 



promise (Hsu & Putman 2008). 



Expanding upon our one-object proof-of-concept analysis of Leo T, we also intend to 
do a similar search for dwarf-like clouds in the GALFA-Hl survey data. Any compact clouds 
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lacking associated infrared flux would be very interesting candidates for optical follow up as 
gas-rich local-group dwarf candidates. 

It is also crucial to determine whether the conjecture that motivated this work, that 
lack of infrared emission in clouds is correlated to halo membership, is indeed correct. Now 
that we have at least two LVCs with no discernible dust, we wish to determine whether 
they indeed lie significantly beyond the disk of our Galaxy {\z\ > 1 kpc). Observations are 
underway to examine stars in the lower halo along the line of sight to these objects to look 
for absorption features. In addition, we would like to examine other LVCs that do exhibit 
infrared emission, to determine whether they are near the disk of the Galaxy, as we would 
expect from their dust emission. 
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